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Abstract 

Deep inelastic electron-photon scattering is studied using e'^e" data collected by the 
OPAL detector at centre-of-mass energies ^/WZ ~ M^o. The photon structure function 
F2{x, Q"^) is explored in a range of 1.1 to 6.6 GeV^ at lower x values than ever before. 
To probe this kinematic region events are selected with a beam electron scattered 
into one of the OPAL luminosity calorimeters at scattering angles between 27 and 
55 mrad. A measurement is presented of the photon structure function F2{x,Q'^) at 
(g2) = 1.86 GeV^ and 3.76 GeV^ in five logarithmic x bins from 0.0025 to 0.2. 
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The hadronic structure function of the photon has been measured in deep in- 
elastic electron-photon scattering at various e"'"e~ experiments ( and references 
therein). In the previous OPAL analysis F2{x, Q"^) was measured using tagged elec- 
trons with a range of scattering angles from 60 mrad up to 500 mrad to the initial beam 
direction. This letter describes an analysis of data collected in 1993 and 1994 at e"'"e~ 
center-of-mass energies between 89.2 GeV and 93.2 GeV, using the OPAL silicon tung- 
sten (SW) luminometer as the electron tagger with a clear angular acceptance from 
27 mrad to 55 mrad. This corresponds to a lowest value of 1.14 GeV^ in the selected 
sample, close to the lower limit at which perturbative QCD can be expected to work. 
The energy of the LEP beams is higher than at any previous e"'"e~ collider so lower 
values of the scaling variable x can be reached at any given Q^. The results reported 
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here go down to x = 0.0025. This kinematical region is of particular interest because 
theoretical predictions differ significantly from each other [^|-[Tl[| and because in this 
region the proton structure function is observed to start to increase as x decreases 



In the singly-tagged regime, with one scattered electron tagged in the detector, 
the two photon process can be regarded as the deep inelastic scattering of an e^ with 
four-momentum on a quasi-real photon with four-momentum p. The flux of quasi- 



real photons can be calculated using the Equivalent Photon Approximation [T^. The 
cross-section for deep inelastic electron-photon scattering with a hadronic final state X 
is expressed as [^, 



dxdg2 xQ^ 



where = —q^ is the negative four-momentum squared of the virtual photon and 
a is the fine structure constant. The usual dimensionless variables of deep inelastic 
scattering, x and y, are defined as 

y^^, (2) 
2 p ■ q p ■ k 

In the kinematic regime studied here {y^ ^ 1) the contribution of the term proportional 
to the longitudinal structure function F^{x, Q^) is small and is therefore neglected. 
The structure function formalism of deep inelastic scattering implies that the virtual 
photon behaves as a pointlike probe. However, at very low it can also show a 
hadronic structure. For this analysis, such a hadronic contribution is neglected. 



The scattered electron is detected with the silicon tungsten calorimeters (SW) |Tf 
that are placed around the beam pipe at a distance of 2.4 m in 2; from the interaction 
point on both sides of the OPAL detector, covering polar angles 6 from 25 to 59 mradQ. 
They are cylindrical sampling calorimeters consisting of 19 layers of silicon detectors 
interleaved with 18 layers of tungsten, equivalent to a total of 22 radiation lengths. 



For electromagnetic showers an energy resolution of 24%/a/-E' [GeV] and a resolution 
in polar angle of 0.06 mrad is achieved. The hadronic final state X is measured with 
SW, the forward detectors FD, the OPAL electromagnetic calorimeter and the OPAL 
tracking system which consists of a silicon microvertex detector, a drift chamber vertex 



detector, a jet chamber and z-chambers [16 



The measurement of F2{x, Q^) involves the determination of and x that can be 
obtained from the four-vectors of the tagged electron and the hadronic final state: 



Q"^ ^ 2 E\j Etag (1 — cos6'tag) (neglecting the electron mass) (3) 



right-handed coordinate system is used. The x-axis points towards the centre of the LEP ring, 
the y-axis upwards and the z-axis in the direction of the electron beam. The polar angle 6 and the 
azimuthal angle (j) are defined with respect to the z-axis and x-axis, respectively. 



4 



i^tag and 6'tag are the energy and polar angle of the observed electron, i?b is the beam 
energy, and W the invariant mass of the hadronic final state. In addition to the tag 
requirement, an antitag condition is applied to ensure that the virtuality of the quasi- 
real photon is small and requirements are imposed on the hadronic final state to 
reject residual background. 



1. A tagged electron is required, identified as a cluster in the silicon tungsten 
calorimeter with energy 0.775 £'b < -Etag < 1-2 ^'b and a polar angle 27 < 6'tag < 
55 mrad with respect to the beam axis, defining a lower limit of Q^in" ^-^^ GeV^ 
and an upper limit of <5max= 6.57 GeV^ 

(tag requirement). 

2. The energy E'at of the most energetic cluster in the hemisphere opposite to the 
one containing the tagged electron is restricted to i?at < 0.25 i?b 

(antitag requirement). 

3. The visible invariant mass W^vis of the hadronic system is required to be in the 
range 2.5 GeV < W^s < 40 GeV. W^s is defined by all tracks of charged parti- 
cles and all calorimeter clusters which are not associated with tracks, including 
clusters from FD and SW but excluding the tagged electron. The same quality 
criteria are applied to all calorimeter clusters and charged tracks as in the previ- 
ous analysis [|l|. The masses of all particles in the hadronic system are assumed 
to be equal to the pion mass. 

4. Nch > 2, where N^h is the number of charged particle tracks originating from the 
hadronic final state. 

5. The transverse momentum component pJi"*' of the event parallel to the tag plane 
has to be lower than 3 GeV. The tag plane is defined by the momentum vectors of 
the incoming beam electron and the tagged electron. The transverse momentum 
of the event is the vector sum of the momenta perpendicular to the beam axis of 
all measured particles, including the tagged electron. 

6. The transverse momentum component of the hadronic system perpendicular to 
the tag plane has to be lower than 3 GeV. 



These requirements select 7112 events, corresponding to an integrated luminosity of 
70.8 ± 0.2 pb-\ with {Q^)= 2.8 GeV^. Using sets of independent trig gers, the trigger 
efficiency was evaluated to be (98±2)%. 

The dominant background comes from two-photon events with a lepton pair in the 



final state. These events are simulated with the Vermaseren Monte Carlo generator |T7 
Background from other sources such as other QED processes with four fermions in 
the final state or the process Z° hadrons is found to be below 0.5%. The total 
estimated background contamination in the selected sample is (2.5 ± 0.2)%. Deep 
inelastic electron-photon scattering is simulated with the generators HERWIG 5.18d 
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|T§, PYTHIA 5.722 |TP[ and F2GEN |2g. F2GEN includes two final state models, 



■'pointlike" and "peripheral", to describe the hard and soft limit of the process (see 
Ref. [|I| for details). All Monte Carlo events are passed through a detailed detector 



simulations program and the same reconstruction and analysis chain as the real 
data events. 

In Ref. [0] several distributions of measured final state hadronic quantities were 
studied in the range 4.6 GeV^ < < 30GeV^, and significant differences were 
observed, both between data and the Monte Carlo models, and between the different 
Monte Carlo models. These distributions have been re-examined in the different kine- 
matic range of this analysis and a similar pattern of disagreements is seen, especially 
in the energy flow distributions (Figure 0), and in the summed transverse energy per- 
pendicular to the tag plane. These differences are found to be most prominent for 
Xvis < 0.05, and small at higher Xvis- (The peripheral final state model of F2GEN has a 
similar behaviour to PYTHIA and is not shown.) Such uncertainties in the modelling 
of the hadronic final state give rise to uncertainties in the unfolded F2{x, Q^) [Q. They 
are included in the systematic error using the same set of Monte Carlo models as in 
Ref. [^, except for the F2GEN models. For F2GEN, a weighted mixture of the two 
final state models available in the generator is used throughout the analysis. This mix- 
ture has been optimised to improve the final state description by a fit to the hadronic 
energy flow of the data in the lowest x bin. To account for the uncertainty of this 
'ad hoc' procedure a variation of the fit by ± three times its error is included in the 
evaluation of the systematic error. 



The program RUN by Blobel p2i is used to unfold the structure function from 



the measured Xvis distribution. To resolve very low x values the unfolding is performed 
on a logarithmic x scale. 

The ability of the unfolding to recover the underlying structure function F2 of the 
data is tested by unfolding the known structure function of Monte Carlo samples, the 
"mock data" , instead of measured data. These mock data samples are then unfolded 
with other Monte Carlo sample. The unfolded F^ functions are compared to the origi- 
nal structure functions used in the generation of the mock data sample. Figure ^(a) and 
(b) show unfoldings of a mock data sample generated with the HERWIG generator and 
the GRV-LO parton density functions using unfolding Monte Carlo samples from 
different generators and different input structure functions. The samples are divided 
into two bins with (a) < 2.5 GeV^ and (b) > 2.5 GeV^, respectively. The er- 
ror bars include statistical errors only. The solid line is the GRV-LO structure function 
for the average of the mock data sample. Figure |^(c) and (d) show similar plots 
for a mock data sample generated with HERWIG and the DC |^ structure function. 
Here, the solid line is the DC structure function. 

The deviations observed represent the systematic impact on the unfolded result of 
different parton density functions and modelling of the hadronic final state. The DG 
structure function is used as a test for an F2 function that vanishes for x ^ 0. The DG 
parton density functions have been evolved by their authors from =4.0 GeV^ 
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and are not supposed to be valid at lower Q^, but they are suitable for this purely 
technical purpose. Figure 0(a) and (c) demonstrate that despite the systematic errors 
originating from different Monte Carlo models a structure function falling for x ^ 
could be measured and distinguished from the GRV-LO prediction in the lower 
bin. In the higher bin the systematic errors are larger: The spread of the unfolding 
results shown in Figure 0(d) does not allow the exclusion of the function described 
by GRV-LO. 

To determine the central values of the measured a "reference" unfolding is de- 
fined. It is based on a HERWIG sample generated using the GRV-LO parametrisation, 
chosen for consistency with previous OPAL results The event selection as de- 
scribed above is applied. The data are divided into two Q'^ bins with Q'^ < 2.5 GeV^ 
and > 2.5 GeV^ which are unfolded separately. The two bins contain approximately 
equal numbers of events. The x binning is chosen to keep correlations low between the 
unfolded values in the different x bins. 



The Monte Carlo generators HERWIG and F2GEN predict mean values of (P^) = 
0.03 — 0.08 GeV^ for the virtuality of the quasi-real photon, depending on the model 
and the structure function used. Several theoretical predictions exist for how F2 should 
behave as a function of [p^ -|25[| . An estimate for the effect of the non-zero virtuality 

yields an increase of F2 by roughly 10%, based on the P^ dependent structure 



function parameterisation of Schuler and Sjostrand |^ and the reference Monte Carlo 
sample. As the distribution of P^ in the data and the correct theoretical prescription 
are not known, no correction is applied to the results. 



[GeV^] 


[GeV2] 


a 

-logio(x) 


: range 

X 


X (centre of 
logio(x) bin) 


F^/a 


1.1 - 2.5 


1.86 


2.6 - 2.2 
2.2 - 1.7 
1.7- 1.4 
1.4 - 1.0 


0.0025 - 0.0063 
0.0063 - 0.020 
0.020 - 0.040 
0.040 - 0.100 


0.004 
0.011 
0.028 
0.063 


0.27 ±0.03^°:°^ 
0.22 ±0.02+°:°2 

0.20 ±0.02+°:°^ 
0.23 ±0.021°:°^ 


2.5-6.6 


3.76 


2.2 - 1.7 
1.7 - 1.4 
1.4 - 1.0 
1.0 - 0.7 


0.0063 - 0.020 
0.020 - 0.040 
0.040 - 0.100 
0.100 - 0.200 


0.011 
0.028 
0.063 
0.141 


0.35 ±0.03^[!:J^^ 
0.29 ±0.03^°:°^ 
0.32 ±0.02^°:°^ 
0.32 ±0.03^|];J^^ 



Table 1: Results for F2 as a function of x in bins of Q^. The first errors is statistical 
and the second systematic. The systematic errors do not contain the systematic effect 
caused by P^ being different from zero. 



The unfolded F2 /a for the data is listed in Table |^, and shown in Figure The 
value of F2 /a is given at the centre of each bin in logio(a;). The error bars show 
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both the statistical error alone and the quadratic sum of statistical and systematic 
errors. Figure ^ also shows the -Fg^/a calculated from the GRV-LO and the SaS- 
ID [[7| leading order parton density parametrisations and the higher order GRV-HO 
parametrisation, evaluated at the corresponding (Q^) values. The parton density 
parametrisations given in Ref. are not shown because they are not supposed to 

be valid at these values, similarly to the DG parametrisation. In addition results 
from PLUTO and TPC/27 [|] for similar values are shown for comparison. 

The central values and statistical errors of the F2 measurements are determined us- 
ing the reference unfolding. The statistical errors include, in addition to the statistical 
error of the measured data, also the error due to the limited number of Monte Carlo 
events (1.6 times the number of data events) which is estimated to be l/Vl-6 times the 
statistical error of the measured data. The estimation of the systematic errors includes 
the following four components: 

• The selection requirements are varied in order to change the signal and back- 
ground event composition and to take into account possible uncertainties in the 
simulation of variables which are used for the event selection. The size of the vari- 
ations reflects the resolution of the measured variable and fulfils the requirement 
that the mean (Q^) of the sample is not shifted significantly by the variation. 

• An altered set of quality criteria for calorimeter clusters and tracks is used to 
determine systematic errors resulting from imperfections in the simulation of the 
detector acceptance and calibration for tracks and calorimeter clusters. 

• The unfolding is performed using the HERWIG generator and the standard selec- 
tion, but replacing the GRV-LO parton density functions with SaS-lD in order 
to study the uncertainty due to the structure functions assumed in the Monte 
Carlo samples. 

• The unfolding is performed with the standard selection requirements using 
PYTHIA and F2GEN in order to study the effect due to a different modelling of 
the hadronic final state in the different Monte Carlo programs. 

For each of the four systematic studies, the maximum deviations (above and below) 
of the various unfolding results from the result of the reference unfolding are taken as 
systematic errors. The total systematic error assigned to the results in Table |I] is the 
quadratic sum of these four contributions. 

In summary, OPAL data recorded in 1993 and 1994 have been used to measure the 
photon structure function F^(x, Q^) at low x using events with an electron tagged in 
the silicon tungsten calorimeters. -F^(x,Q^) has been unfolded as a function of x in 
two bins of Q^, with (Q^) = 1.86 and 3.76 GeV^. The unfolding has been performed 
on a logarithmic x scale to resolve the lowest accessible x values. The data have been 
unfolded down to a minimum x of 0.0025, lower than measured previously |l|-^. This 
is the region where the proton structure function as measured at HERA starts to 
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rise 



. The extracted F2 result is consistent with the other OPAL F2 measurements 



Within errors the results agree with the lowest value published by PLUTO 0. 
Compared to the TPC/27 measurement our results tend to be higher. Our result 
is consistent with a flat -F^(x) in both ranges within the errors though it does not 
exclude a small rise with decreasing x. The unfolded result is consistent in shape with 
the GRV-LO and SaS-lD parameterisations for the corresponding values. However, 
the measured is higher than the GRV-LO and SaS-lD predictions. The GRV-HO 
prediction follows the data more closely and is in good agreement in the lower bin. 
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Figure 1: The hadronic energy flow per event as a function of pseudorapidity i] for 
tfie data and various Monte Carlo samples. In (a)-(f) the values of the bin limits in 
Xvis have been derived from the bin limits in Xtruc of the three lowest x bins in Figure 
^(a). The bins in 6'tag correspond approximately to the bins of the unfolding. The 
errors shown are statistical only. The vertical lines show the acceptance regions of the 
OPAL detector components, CD = Central Detector, FR = Forward Region and BP 
= Beam Pipe. (FR, BP not marked on the tag side). 
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Figure 2: Unfolding tests: In (a) and (b) mock data samples of HERWIG events 
with the GRV-LO parton density functions were unfolded using one of the four listed 
unfolding Monte Carlo samples. The solid curves show the GRV-LO F'^ix) for — 
1.86 GeV^ and = 3.76 GcV^. In (c) and (d) a similar exercise was performed 
with mock data samples from HERWIG with the DG parton density functions. The 
solid curves show the corresponding DG values for F2{x) at = 1.86 GeV^ and 
Q'^ — 3.76 GeV^. The error bars are statistical only. The symbols are slightly shifted 
in X to avoid overlap. 
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Figure 3: The full circles show our result for ^2(0;) at (Q^) = 1.86 GeV^ (a) and 
(Q^) = 3.76 GeV^ (b). The total error and the statistical contribution are shown for 
each point. The tick marks at the top of the figure show the bin hmits in x for both 
ranges. These points are placed in the middle of the bin. The curves indicate the 
GRV-HO (dotted), GRV-LO (solid), and SaS-lD (dashed) predictions for at the 
corresponding Q^. The open symbols show results from other experiments at similar 
Q^, with the total errors indicated only. These points are placed at the centre of the 
bins on a linear x scale. 
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